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ABSTRACT 

We consider implications of a possible presence of a Thomson-thick, low-temperature, plasma 
cloud surrounding the compact object in the binary system Cyg X-3. The presence of such a 
cloud was earlier inferred from the energy-independent orbital modulation of the X-ray flux 
and the lack of high frequencies in its power spectra. Here, we study the effect of Compton 
scattering by the cloud on the X-ray energy and power spectra, concentrating on the hard 
spectral state. The process reduces the energy of the high-energy break/cut-off' in the energy 
spectra, which allows us to determine the Thomson optical depth. This, together with the 
observed cut-off' in the power spectrum, determines the size of the plasma to be ~ 2 x 10^ 
cm. At this size, the cloud will be in thermal equilibrium in the photon field of the X-ray 
source, which yields the cloud temperature of ^ 3 keV, which refines the determination of 
the Thomson optical depth to ~ 7. At these parameters, thermal bremsstrahlung emission of 
the cloud becomes important as well. The physical origin of the cloud is likely to be collision 
of the very strong stellar wind of the companion Wolf-Rayet star with a small accretion disc 
formed by the wind accretion. 

Key words: accretion, accretion discs - binaries: general - radio continuum: stars - stars: 
individual: Cyg X-3 - X-rays: binaries - X-rays: stars. 



1 INTRODUCTION 

Cyg X-3 is a high-mass binary with a Wolf-Rayet companion 
jvan Kerkwiik et with an unusually short orbital pe- 

riod of 4. 8 h, located a t a distance d — 9 kpc in the Galac- 
tic plane jPickevI 1 19831 : IPredehl et all I2OO0I) . Due to the lack 
of reliable mass functions and determination of the inclina- 
tion, it remains uncertain w hether its compact object is a black 
hole or a neutron star (see I Vilhu et alj [2OO9I for a recent dis- 
cussion). However, the presence of a black hole is strongly 
favoured by considering the X-ray spectra, the radio em ission, and 
the b olometric luminosity 1 Szostek & ZdziarskH I2OO8I hereafter 
SZ08: ISzostek. Zdziarski & Mc Collough 2008,, h ereafte r SZM08; 
iHialmarsdotter et alj l2008l : iMalmarsdotter et afl l2009l hereafter 
Hj09). 

Cyg X-3 is a persistent X-ray source. Its X-ray spectra have 
been classified into 5 states by SZM08, who have also quantified 
their correlations with the radio states. Fig. [T] compares the char- 
acteristic spectra of a typical black-hole binary, XTE J 1550-564 
jZdziarski & Gierlinskill2004l) . with the states of Cyg X-3. We see 
a remarkable overall similarity between the spectra of the two ob- 
jects, apart from the very strong absorption in Cyg X-3. Its spec- 
trum 5 is similar to the ultrasoft one of XTE J1550-564. Then, the 



* E-mail: aaz@camk.edu.pl (AAZ), rmisra@iucaa.ernet.in (RM) 



spectra 3 and 4 are similar to the very high and high, respectively, 
states of XTEJ1550-564. Finally, the hardest spectra of Cyg X-3, 1 
and 2, can be identified with the canonical low/hard state of black- 
hole binaries. Here, however, the high-energy cut-off/break occurs 
in Cyg X-3 at a lower energy, ~ 20 keV, than in other black-hole 
binaries, where it is at ~ 100-200 keV, as seen in Fig.[T] 

Initially, X-ray spectr a of Cyg X-3 were mo d elled rather phe- 
nomenologically, e.g., by IWhite & Holtl ( Il982h . iNakamura et al.l 
( ll993h , lRaieev et al.l ( fT994h . including, e.g„ power-law, blackbody 
and/or bremsstrahlung components. M ore recent ly, physically mo- 
tivated models have been studied by IVilhu et al. ( 2003), SZ08, 
IHialmarsdotter et al] ( l2008h and Hj09. When the hard-state spectra 
were fitte d by thermal Comptonization, the dominant process in the 
hard state jZdziarski & Gierlinskil2004l) . the obtained electron tem- 
perature was kT\i ~ 4-8 keV, much lower than kT\i ~ 50- 100 keV 
in the usual hard state of black-hole binaries (e.^.,,Gierlinski et al.l 
ll997l : IZdziarski et al.lll998l : IWardzinski et al.ll2002l) . 

An important, but not included yet in the spectral stud- 
ies, physical process in Cyg X-3 may be Compton scattering 
by a relatively cold, ionized, plasma with a Thomson optical 
depth > 1, surrounding the X-ray source in Cyg X-3. The 
presence of a highly ionized Thomson-thick absorber is indi- 
cated by a deep ~9 keV Fe K edge seen in the X-ray spec- 
tra (see, e.g., SZ08). The presence of a Thomson-thick plasma 
in Cyg X-3 was postulated by, e.g., iDavidsen & Ostrike^ ( Il974h 
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Figure 1. (a) The canonical spectral states of black-hole binaries as exhibited by XTE J 1550-564 jzdziarski & Gierlihskill2004 . (b) The main spectral states 
of Cyg X-3 as classified by SZM08. Identical colours on both panels denote the corresponding states, except that both the blue (1) and cyan (2) spectra of Cyg 
X-3 c orrespond to the low/hard state. The spectra are modelled (solid curves) by hybrid (composite thermal and nonthermal) emission jZdziarski & Gierliiislj 
|2004 SZM08); the dashed curves in (a) show the spectra before absorption. 



and iHertz. Joss & Rappaport 1 1978 ) to account for the X-ray or- 
bital modulatimrAilojBe^^ derKlisl ( ll994h proposed the 
presence of such plasma based on the power spectra of Cyg X- 
3. That study, however, was hampered by an unresolved (at that 
time) instrumental effect, which prevented iBerger & van der 

Mi 

i ll 9941) from performing a more detailed study of the effect of 
scattering on the power spectra. The lack or weakness of X-ray 
variability at frequencies of / > 0. 1 Hz in Cyg X-3 was con- 
firmed by .Axe lsson, Larsson & Hialmarsdotter l 2009h. We note 
that neitherloavidseri & Ostrikeil jl974l) . iHertz et al-l d 19781) nor 
iBerger & van der KUsI l ll994h discussed effects of the scattering 
cloud on the X-ray spectra of the source. This elfect was mentioned 
by Hj09, but not taken into account in their fits due to the lack of a 
readily available spectral model. 

The location of the T ho mson-thick plasma re mains unclear. 
IVan Kerkwiik et al.l ( 1 19921) and lvan KerkwiikI ( Il993h estimated the 
Thomson optical depth of the stellar wind from the Wolf-Rayet 
companion measured from the compact object towards the observer 
to be ~ 3. On the other hand, SZ08 found it to be < 1. G enerally, 
as discussed by, e.g., SZ08 or Ivan Kerkwiik et all ( Il996l) . the esti- 
mates of the mass-loss rate and the wind clumping in Cyg X-3 are 
uncertain, making these estimates relati vely uncertain. Th e wind 
may also form a circumbinary envelope dVilhu et al]|2007b . which 
may be optically thick. On the other hand, outflows from an accre- 
tion di^c_aroundthe compa ct object may be optically thick as well, 
e.g.. |Poutanen et al] ( l2007t) . However, such outflows are usually ax- 
ially symmetric around the compact object, which would not result 
in the observed orbital flux modulation. Another possibility for the 
physical identification of the scattering cloud is a bulge formed by 
the stellar wind colliding with the accretion disc. 

In this paper, we study in detail the effect of the presence of 
such plasma surrounding the X-ray source, on both the X-ray spec- 
tra (Section |2} and on the power spectra (Section O in the hard 



state. In Section |4| we present our final model reproducing the en- 
ergy and power spectra of Cyg X-3 and satisfying thermal equilib- 
rium in the scattering cloud. In Section [5] we discuss our results 
and possible application of our model to the soft states of Cyg X-3. 



2 COMPTON DOWN-SCATTERING OF X-RAYS 

We consider the average hard-state spectra 1 and 2 of SZM08 
from the RXTE Proportional Counter Array (PCA) and the High 
Energy X-Ray Transient Experiment (for the log of the observa- 
tions, see Hj09; note that those authors used the term 'interme- 
diate' for the spectrum 2). The advantage of choosing the hard 
state is that it has usually a standard form with a pronounced high- 
energy cut-off, and it is well modelled by thermal Comptonization 
in a hot plasma (e.g., Gierliiiski et al.*'l997'; Zdziars ki et al. I I 19981 : 
IWardzinski et al. 2002; Zdziarski & Gierlinski 2004). On the other 
hand, high-ener gy tails in soft states exhibit a variety of spectral 
shapes, see, e.g., IZdziarski & Gierliiiskil ( l2004l) . 

We use the therm al Comptonization model of 
IPoutanen & Svenssorj jl99^ in spherical geometry (with the 
geometry parameter = 0), compps in xspec (Arnaud 1996). As 
the parameters of the model we use the electron temperature, Tj,, 
and the Compton parameter, = 4Thfcrh/meC^, where th is the 
Thomson optical depth of the hot plasma, and is the electron 
mass (which definition is suitable for Th ^ 1). In addition, we 
include Compton reflection (Magdzi arz & Zdz iarski 1995) from 
a highly ionized medium with an associated Fe K line, assuming 
a reflector solid angle of 2n at an inclination of 60° and the line 
equivalent width of 100 eV. (Although we find below that the 
reflection features in the intrinsic spectra are strongly smeared out 
in the scattered spectra, we keep them for physical consistency.) 
We assume seed photons Comptonized by the hot plasma to have 
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19881: IZdziarski, Johnson & Magdziar^ 11996 

Zdziarski. Poutane n & Johnson 200G; Zdzi arski et alj l2003ll 
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Figure 2. Effect of Compton down-scattering on X-ray spectra in Cyg X- 
3. Tlie error bars show tlie average liard-state spectra, groups 1 and 2 of 
SZM08, the dotted curves show the incident spectra, and the solid curves 
show the spectra after transmission through a sphere of cold ionized elec- 
trons with the optical depth of tq. The incident spectra are from thermal 
Comptonization in a hot plasma with an electron temperature, kT[i, and the 
Compton parameter, v'h, and including Compton reflection corresponding to 
a solid angle of an ionized reflector of 2n with an associated Fe K line. The 
crosses show the transmitted spectra after absorption by an external mat- 
ter with a column density, Nh, and including an externally produced Fe K 
line, (a) The average spectrum 1 compai'ed to the down-scattered thermal 
Compton spectrum with kTi, = 17 keV, Vh = 1-4, tq = 4.5, Nu = 9 X 10^^ 
cm"^. (b) The spectrum 2 modelled by kTii = 25 keV, yi, = 0.9, To = 5.0, 
A^H = 3.5 X 10^^ cm-2. 



a disc blackbody spectrum with the maximum temperature of 0.3 
keV. 

As shown e.g., by SZ08, intrinsic X-ray absorption in Cyg 
X-3 is very complex, caused by the stellar wind with at least two 
phases. Here, we do not attempt to reproduce this complexity as it is 
not related to the main goal of this paper. Instead, we simply model 
absorption (of the photons leaving the scattering cloud) as due to a 
neutral medium. Also, we do not physically model the observed Fe 
K line, produced apparently outside of the scattering cloud. Instead, 
we simply add a Gaussian line reproducing the observed line to our 
final spectra. 

Since hard-state spectra of black-hole binaries extend well 
above 100 keV, we need to accurately take into account the 
Klein-Nishina effects in the postulated scat tering medium. For that 

f urpo se, we us e the Monte Carlo code of iGorecki & Wilczewskil 
1984) and loierl inski ( 2000). The c ode has been extensively tested 
against other codes and solutions jwhite, Lightman & Zdziarsla 



We use a spherical geometry with uniform electron density, n^, 
and the source of photons at the centre. We note that the scattering 
medium in Cyg X-3 is likely to be spherically asymmetric if it is 
formed by accretion, and it needs to be such in order to account 
for the observed orbital modulation (Section [TJ. However, here 
we attempt to reproduce only the properties averaged over the 
orbit, in particular, the average X-ray spectra of SZM08, for which 
purpose the assumed spherical symmetry is fully sufficient. The 
input photons have the X-ray spectrum as described above. In this 
section, we consider only electrons at a temperature low enough 
for the effect of Compton up-scattering to be negligible. This 
happens when the Compton parameter of the scattering medium, 
yo = 4(kTo/m^c^)Tg (which form of this definition is valid for 
To s> I), is ^ I, where To is the temperature, tq = n^o-jR is the 
Thomson optical depth, R is the radius, all of the scattering cloud, 
and (Tj is the Thomson cross section. In this section, we assume 
kTo = 0.2 keV, which satisfies the above condition for the values 
of To we find here. 

The models we consider are relatively complex. In particular, 
they involve a convolution of intrinsic spectra produced by ther- 
mal Comptonization in a hot plasma with the effect of Compton 
scattering in the, much colder, surrounding cloud. In Section|4]be- 
low, we also find varying of the cloud temperature to be impor- 
tant and Comptonization of thermal bremsstrahlung to be a ma- 
jor process. At present, convolution models describing Compton 
scattering in a plasma of arbitrary Thomson optical depth are not 
available in xspec. Given the complexity of this problem, we defer 
obtaining Green's functions for Comptonization to a future work. 
Thus, we are unable to formally fit X-ray count spectra. Instead, we 
compare our models to the X-ray spectra of SZM08, which were 
de-convolved assuming hybrid-plasma models. This will only en- 
hance some discrepancies of the present models with the spectra of 
SZM08 obtained using different models, but it would not lead to a 
false agreement. 

We have found that we can reproduce the average hard-state 
spectra with tq - 5, fcTh ~ 20 keV and Vh ~ 1, see Fig. |2] 
for details. The position of the peak of the down-scattered spec- 
trum depends only weakly on kTt,, as it is roughly at m^c^lT^ 
dSunvaev & Titarchukll98a) : however, the slope of the high-energy 
tail of the spectrum depends sensitively on it. The hardness of the 
spectrum below the peak depends on the Compton parameter, 
therefore we obtain a lower value of it for the softer spectrum. As 
seen in Fig. [2l we can reproduce rather well the high-energy parts 
of the spectra, >10 keV. On the other hand, as discussed above, 
accurate modelling of the part of the spectrum at <10 keV would 
require taking into account complicated physical processes in the 
inhomogeneous stellar wind, which we do not attempt here. 



3 THE EFFECT OF COMPTON SCATTERING ON 
TIMING PROPERTIES 



iLightman. Lamb & Rvbickil ( Il98lh calculated the distribution of 
the number of scatterings experienced by a photon emitted at the 
centre of an optically- thick sphere of c old electrons, given by their 
equation (18). K vlafis & Klimisl ( ll987h used that result to calculate 
the distribution of the time spent by a photon inside the sphere, 
AN I At, with certain additional approximations. In particular, they 
assumed that a photon scattering n times has covered the Thom- 
son optical depth of exactly r = n along its path, i.e., they re- 
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placed directly « by r in equation (18) of lLightman et alj ( Il98lh . 
On the other hand, the distribution of the number of scatterings 
af ter a given elap s ed tim e is Poissonian, see, e.g., equation (6) 
in llllarionov et al.l ( Il979l) . Thus, we can sum over n for a given 
elapsed time, improving the accuracy with respect to the results of 
lKvlafis&Klimisl ( ll987h . In particular, we find. 
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for t>0.5tq. Here, t is both the optical depth that a photon has 
travelled and the elapsed time in units of the Thomson time, r = 
ctT(,/R, where R is the sphere radius. 

We have compared these results against those from our Monte 
Carlo code in the spherical geometry. We have considered photon 
energies in the Thomson regime since power spectra measured by 
RXTE PCA correspond mainly to it. We note here that there is a 
difference between the distribution of the photon arrival time at the 
sphere surface and that measured by a remote observer. Apart from 
a constant time shift, the latter is equal to the distribution measured 
at a plane tangent to the sphere. In units of the Thomson time, the 
difference is At = Tq(1 - cos a), where a is the angle at which a 
photon leaves the sphere with respect to the radial direction. Thus, 
At = for photons leaving the sphere radially, and Tq for photons 
tangent to the sphere surface. 

Usi ng the Monte Carlo co de, we have found that the results 
based on lLightman et alj il98lh approximate the actual distribution 
only for a very large optical depth of the sphere, Tq > 20 or so. Even 
including summation over the Poissonian distribution gives still a 
poor agreement with the Monte Carlo results at Tq = 5, 7.5, as we 
show in Fig. [3] 

We have thus obtained a set of fitting formulae for dN/ dr for 
values of Tq of interest for this work. If normalized to unity, they 
also give the Green's functions, G(t) = dN/dr, 



G(T): 



exp[-^-(^)'] + 5(T-To)e- 



T < To; 
r > To, 



(2) 



Here, Ti, T2 are fitted, the term with the delta function accounts 
for the un-scattered photons, and r follows from the normalization 
condition of J^°° G(T)dT = 1. For To = 5 and 7.5, t, ^ 9.513, 19.87, 
T2 7.231, 14.61, r ^ 0.2813, 0.1299, respectively. We see in 
Fig. [U that equation ^ provides excellent fits to the Monte Carlo 
results. The actual distributions are significantly flatter than those 
corresponding to the results for tq » 1 . In particular, equation (T} 
for To = 5 yields ti ==8.11 (and ti = Sr^/n- =^ 7.60 if the summa- 
tion over the Poissonian distribution is neglected). 

If the intrinsic source variability is 5(t), the observed light 
curve, F(t), is a convolution, 



F(t) 



f 



5(T-T')G(T')dT'. 



(3) 



The Fourier transform of F{t), T , is then the product of the indi- 
vidual transforms, i.e.. 



T(0) = S{<f>)g{<f>), g(<f>) : 



£ 



G(T)e 
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where (p = /R/tqc is the dimensionless frequency in the Thom- 
son units. If the original signal is sinusoidal, its amplitude is re- 
ducecQ by \g{(p)\. In general, damping due to scattering of an in- 
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Figure 3. The distribution of the time spent by a soft photon in a spheri- 
cal cloud of cold electrons measured from the cloud centre in units of the 
Thomson time for (a) To = 5, (b) tq = 7.5. The black dots give Monte Carlo 
resuhs and the red soli d curve gives its fit, equ ation (2). The dotted curves 
gives the distribution of lLightman et alj il98ll) corrected for the Poissonian 
distribution of the time corresponding to a given number of scatterings. 



trinsic variability can be calculated either using the convolution 
of the light curves, equation {3), or by integrating the resulting 
power spectrum, equation l|4j. The intrinsic power spectrum at 
is damped by \g{(p)f, which quantity is shown in Fig. |4| We note 
that our Monte Carlo results show that the Fourier transform ob- 
ta ined bviKvlafis & PhinnevI ( [l98ft) . which is based on the results 
of lLightman et alj ^ 1981 ), becomes accurate only at To ^ 20. 

Based on our Monte Carlo results, we have found that the 
characteristic dimensionless frequency, (p^, which we define as 
that at which scattering reduces \g((/))f by a factor of e, is (p^ - 
1/(2tq), implying a characteristic dimensional frequency of fc - 
c/(2Rto). The characteristic frequency is 4 times lower than that 
obtained using the inverse of the (commonly used) average num- 
ber of scattering in a sphe re with a central source, Tq/2 (e.g., 
ISunvaev & Titai^chukl Il980l) . Thus, our quantitative calculations 
demonstrate that Compton scattering is significantly more effec- 
tive in damping variability than commonly assumed, affecting fre- 
quencies ~ 4 times lower than thought before. This is due to the 
relatively broad distribution of the number of photons per InT 
(= TdN/dr). 

We then compare our results with the observed power spectra 
of Cyg X-3. We have computed them for the hard-state RXTE ob- 
servations used for the average hard-state spectrum 1 of SZM08, 
and the resulting average spectrum is shown in Fig.jS] We see a dis- 



We note that we need to transform the full Green's function, including transfo rm of the scattered photons, which was done in iKvlafis & PhmnevI 

the un-scattered part, rather than to add e"'"" to the absolute value of the (l98i). 
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Figure 4. The squared attenuation coefficient as a function of frequency in 
the Thomson units {</> = /R/tqc). The solid, dashed and dotted curves are 
for To = 5, 7.5, 10, respectively. At high frequencies, \G(<f))\ —> e"^". 



10-^ 




Figure 5. The average power spectrum per In/ of Cyg X-3 from the 
RXTE/PCA observations included in the average spectrum 1 (crosses). The 
solid curve shows a fit of for tq = 7.5 to the data at / > 0.04 Hz. 



tinct cut-olf at a frequency of =^ 1 Hz. Using our results above, 
we can then deduce the size of the scattering cloud as 

r 1.5x10'° cm 

(5) 



R 



2ro/c ro(/JlHz) 



For ~ 1 Hz and Tq = 5,R ~ 3x10'' cm, which is <k the separation 
between the components in Cyg X-3, a 3 x 1O"(A^/3OM0)''\ 
where M = Mx+M2, and Mx and M2 are the masses of the compact 
object and the donor, respectively. 



4 THERMAL EQUILIBRIUM OF THE SCATTERING 
CLOUD 

An important consequence of the relatively small size of the scat- 
tering cloud derived above is that its cooling by the compan- 
ion star will be ineffective. The stellar luminosity is estimated as 
L, - (1-30) X 10^^** erg s"' (see, e.g., a discussion in SZ08). Given 
R/a 10"^, the cloud receives only ~ 10"*L», which is <k the X- 
ray luminosity of Lx ~ 10'* erg s"' regardless of the exact value of 
L,. Thus, the scattering cloud will be effectively heated by photons 
of the X-ray source and the surrounding accretion disc, and will 
achieve an equilibrium temperature, kTo, with radiative gains bal- 
anced by radiative losses. The amplification factor. A, of the Comp- 
tonizing cloud can generally be written as. 



where Lj, is the luminosity of the central hot-plasma source (which 
is the only source of energy, and thus = Lx), and La, is the lumi- 
nosity of this source after passing through the scattering cloud. For 
our parameters, A < 1 usually (though not necessarily), depending 
on both the parameters of the cloud (Tq, Tq) and the shape of the 
intrinsic spectrum, determined by Lh and Vh- 

We thus allow a free kTo imposing initially the condition of 
Lch = L\i- However, we find we can satisfy it only very roughly 
while keeping the model spectra fitting the data. The reason is that 
increasing kT^ moves the peak in the down-scattered spectrum up, 
which requires an increase of tq (which moves the peak down) to 
compensate for it, which in turn reduces Lch/Lj,, making it still 
< 1. The obtained temperatures are kTo ~ 3 keV, and the optical 
depths somewhat higher than in the models above, tq ^7. For these 
parameters, we find, howev er, that bremsstrahlung emission (e.g., 
iRvbicki & Lightmad[T979l) of the scattering plasma is a substan- 
tial source of photons at the above kTo and To. The bremsstrahlung 
luminosity, Lj,, of a spherical cloud can be calculated to be. 



2"-n^'^afgqfj,TlRm,c^ j kTo 
1.8xl0"|.(i^f 



R 



1/2 



2 X 10" cm 2 \3keV 



kT, 



1/2 



erg s 



(V) 
(8) 



where fi^ - 2/(1 -1- X) is the mean electron molecular weight, 
X is the mass H relative abundance, cff is the fine-structure con- 
stant, g is the energy-averaged Gaunt factor, and q accounts 
for a density gradient and/or du mpiness in the scattering cloud 
JOwocki. Gavlev & Shavivll2004l) . 



{n,f/{nl) > 1. 



(9) 



For consistency with our treatment of the timing properties, we 
consider a uniform cloud, where q = I. The factor of yUe follows 
from averaging over Z" at a given Tq in a mixture of H and He, 
where Z is the atomic charge. For cosmic composition, X ^ 0.7 and 
fic ^ 1.2. However, Cyg X-3 ha s a Wolf-Rayet companio n, which 
composition is dominated by He jvan Kerkwiik et al.ll996l : see also 
a discussion in SZ08), for which fi^ 2. Thus, bremsstrahlung is, 
for a given To, a factor of 1 + X ^ 1.7 more efficient in a He plasma 
than in a cosmic-composition medium. Note that L\, oc R for a 
given To. We model the bre msstrahlung spectrum using the en ergy- 
dependent Gaunt factor of iKellogg. Baldwin & Kocti ( Il975h . The 
numerical values in equation ^ correspond to our final solution 
for the spectrum 1 below. 

The bremsstrahlung emission is then Comptonized in the scat- 
tering cloud. Unlike Comptonization of the intrinsic central source, 
the seed photons are now uniformly distributed in the cloud (which 
follows from our assumption of the uniform cloud density), which 
we take into account in our Monte Carlo modelling. (We note that 
bremsstrahlung Comptonization could also have been dealt with 
using the Kompaneets equation, though less accurately than by the 
Monte Carl o method given the uniform distribution of the seed 
photons, see lSunvaev & Titarchu k 1980.) Comptonization depends 
(weakly) on bremsstrahlung self-absorption, which provides an ef- 
fective low-energy cut-off for the seed photons. For our parameters, 
the energy at which the optical depth for self-absorption becomes 
unity is 0.3 eV (calculated following Rybicki & Lightm^ 

which value we adopt. The amplification factor of the 
bremsstrahlung emission, Aj,, can be written then as 



■ A(ro,To,rh,jh), 



(6) =Ab(ro,To,LO> 1, 



(10) 
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where Lcb is the luminosity of the Comptonized bremsstrahlung 
emission of the cloud, and the dependence of is weak. 

The Comptonized bremsstrahlung emission by the electrons is 
an important cooling mechanism for the considered problem. The 
energy equilibrium equation can be written as, 

Lh = Lch+Lcb- (11) 

The size of the cloud, R, in equation Q is then determined from 
equation ([5}, which makes Lb oc tq. 

Fitting the observed spectra while at the same time satisfying 
equations (O, dllt is a complex non-linear problem. To facilitate 
understanding of the requirements on the solution, we discuss now 
the effects of varying various parameters. An increase of Tq moves 
the position of the peak in the scattered spectrum to higher energies. 
Also, it increases all Lb, Lcb and Lch. On the other hand, the normal- 
ization of the bremsstrahlung spectrum below its high-energy cut- 
off at ~ kTa decreases as Lq An increase of tq moves the peak of 
the scattered spectrum down, though this saturates when the peak 
reaches ~ SfcTo. Then, it increases Lcb and decreases Lch- Also, it 
steepens the slope of the scattered spectrum beyond its peak. An 
increase of Tb reduces Lcb and makes the high-energy tail harder. 
An increase of Vh makes the slope below the peak harder, which 
holds for both the initial intrinsic spectrum and the scattered one. 
Finally, increasing A'h hardens the low-energy part of the spectrum 
without affecting the energy balance. 

We present the models reproducing the average spectra in Fig. 
(6] which also gives the obtained parameters. We see in Fig.|6]that 
the spectrum 1 can be fitted very well. On the other hand, there is an 
indication of an additional component at low energies in the spec- 
trum 2. Unfortunately, we cannot study it further with the RXTE 
data only, covering only the range > 3 keV. Still, the high energy 
parts are very well reproduced in both cases. The relative contri- 
bution of the Comptonized bremsstrahlung to the total luminos- 
ity is substantial, = 0.3, 0.2 for the spectra 1 and 2, respectively. 
These factors also equal the fractional energy loss of the intrin- 
sic emission in the scattering cloud. Comptonization amplifies the 
bremsstrahlung emission by factors =^ 1.5, 1.3, respectively. We 
also find that including bremsstrahlung makes it possible to obtain 
good fits satisfying exactly the energy equilibrium, unlike the case 
without bremsstrahlung. The energy equilibrium is satisfied with 1 
per cent in the cases shown in Fig. (6) For the obtained parameters, 
the Compton temperature of the scattering cloud is vo ~ 1- Then, 
photons at L > 10 keV are down-scattered, but photons at lower en- 
ergies are efficiently up-scattered by the plasma. 



The obtained temperature of the hot Comptonizing plasma is 
now fcTb = 30-50 keV. Generally, k Tj, ~ 50-100 keV are found i n 
the hard state of black-hole binaries jZdziarski & Gie iiinski'20Ql), 



with lower values found in more luminous states 


Wardzinski etal.l 


l2002l;lYamaoka et alj2006l:lMivakawa et alj2008 


0. In order to see 



if the 50-100 keV range can be extended to lower values, we have 
analyzed the RXTE data for a luminous hard state of the black- 
hole candidate GX 339-4 (from 2002 April 23, the observation ID 
70109-01-04-00), with the bolometric accretion lu minosity of Lx - 
2.4 X 10'^ erg s ' at the assumed distance of 8 kpc jZdziarski et al.l 
|2004|) . and we have found using compps kT^ 33;^2 keV (with 
xl = 82/128 including a systematic error of the data of 2 per cent). 



^ Note that tempera tures as low as ~ 20 k eV were claimed in the hard 
state of GX 339^ bv lMivakawa et al]j2008l) . Howev er, they were obtained 
using the non-relativistic Comptonization model of ISunvaev & TitarchukI 
(1980|), wfiich is not valid for spectra extending to > 50 keV. 




E [keV] 

Figure 6. The observed average spectra (error bars) compared to intrinsic 
thermal Compton spectra (dotted curves) transmitted (solid curves) through 
a Compton-thick cloud at the temperature, Tq, corresponding to thermal 
equilibrium and including the bremsstrahlung emission (dot-dashed curves) 
emitted by the cloud and Comptonized in it (dashed curves). The crosses 
show the sum of the Comptonized intrinsic and bremsstrahlung spectra after 
absoiption by external matter with a column density, A'h and including an 
externally produced Fe K line, (a) The average spectrum 1 compared to the 
down-scattered thermal Compton spectrum with kT^ = 30 keV, = 2.0, 
TO = 7.5, klo = 3.3 keV, A'h = 1.3 x lO^'' cm^^. (b) The spectrum 2 
modelled by kl^ = 50 keV, )'b = 1.2, tq = 7.0, klo = 2.5 keV, A'h = 5 X 
10^^ cm"^. In the thermal equilibrium, the unabsorbed intrinsic bolometric 
flux equals the sum of that flux after Compton scattering in the cloud and 
the flux of Comptonized bremsstrahlung. In both states, the bolometric flux 
is 1.1 X 10"** erg cm"^ s"'. 



On the other hand, the values of Lx implied by our model (at d = 9 
kpc and assuming isotropic emission) in both of the analyzed hard 
states of Cyg X-3 are ^ 1.2 x 10'^ erg s"', relatively similar to that 
in GX 339^. Thus, kTi, ^ 30 keV in Cyg X-3 appears to be an 
entirely plausible value. 

The temporal Green's function for Tq = 7.5 (which is almost 
the same at kTo <k 1 keV as at kTo =^ 3 keV) is approximated 
by equation ^ and it is shown in Fig.[3lb). Its squared transform, 
|^(0)P, is shown in Fig.|4] We have fitted it to the observed power 
spectrum at high frequencies (assuming the intrinsic power oc /"'), 
as shown in Fig. [5] The fit implies R =^ 1.8 x 10' cm, almost the 
same as given by equation l[5). It is somewhat smaller than that 
found in Section |3] due to the lower value of tq - 5 found in the 
models neglecting cloud heating. 
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5 DISCUSSION 

Given the obtained small size of the scattering cloud of ~ 2 x 10' 
cm, it cannot be the stellar wind itself. A likely origin of the cloud 
appears to be collision of the gravitationally focused stellar wind 
with the outer edge of an accretion disc. The resulting cloud or 
bulge will be axially asymmetric with respect to the compact ob- 
ject (as required by the observed orbital modulation of X-rays). The 
presence of such a bul ge was recently found in the high-mass bl ack- 
hole binary Cyg X-1 jPoutanen. Zdziarski & Ibragimovll2008l) . In 
Cyg X-1, the bulge is Thomson thin. However, the orbital separa- 
tion, a, in Cyg X-3 is about 10 times lower than that in Cyg X-1, 
as well as the mass loss rate is at least 10 times higher (e.g., SZ08). 
Thus, the bulge in Cyg X-3 can easily have To ~ 7 and engulf the 
X-ray source. 

The size of the bulge is then similar to that of the disc. Since 
accretion in Cyg X-3 proceeds via stellar win d, its outer radius, 
is mu ch smaller than that of the Roche lobe l lShapiro & Lightmanl 
1 19761) . Using their formalism, ^ 4afo"*(l -I- qf^, where q = 
MilMx, b = Vjei/Vorb, is the relative orbital velocity, v^^j = 
^orb ^wind' ^"'^ "^wini IS the wiud velocity, which is 10** cm s"' 
(e.g., SZ08). Then, ~ R, e.g., =^ 2 x 10'' cm fo r plausible 
values of M = 30Mo and q = 3 jVilhu et alj[2009l) . This also 
yields jShapiro & Lightmanl 1976h the fraction of the mass loss rate 
captured by the compact object in a rough agreement with our value 
of Lx. 

The characteristic ionization parameter of the cloud is ^ ~ 
Lx/intR'^) ~ erg cm s ', which implies that the cloud will be 
almost fully ionized. At the obtained rather high equilibrium kTg ^ 
3 keV, the plasma will also be strongly coUisionally ionized. The 
observed Fe K line and edge are thus formed mostly outside the 
cloud, in the stellar wind (e.g., SZ08). Still, there may be some 
contribution of the scattering cloud to absorption, which may be 
tested with high energy-resolution data, e.g., from Chandra. 

For the sake of simplicity, we have assumed the scattering 
cloud to have a uniform density. In reality, it will have a density 
gradient and no clear edge. An effect of that will be a less gradual 
cut-off of the attenuation function. However, this effect is most im- 
portant for low values of tq. At tq - 7, most of the scattering takes 
place relatively deep in the cloud and the surface effect is minor. 

An important finding of this study is the importance of thermal 
bremsstrahlung, which is otherwise found to be rather unimportant 
in black-hole binaries. In the present case, its importance is related 
to the presence of a dense and compact scattering cloud, apparently 
due to a high density of the stellar wind at a very high mass loss rate 
in this very compact binary. The luminosity of the bremsstrahlung 
turns out to be just at the right level to yield the energy balance of 
the cloud, which is highly remarkable. As we see from equation 
Q, Lb is determined by the cross section of the process, Tq, kTg 
and R, and it could be in principle at any level. Still, the values of 
To and kTo resulting from the spectral fitting and the value of R 
determined by the observed power spectrum result in the value of 
Lb just right to compensate for down-scattering losses of the pri- 
mary emission due to thermal Comptonization in a hot plasma. We 
consider this coincidence a strong argument for the correctness of 
the presented model. We note that an inhomogeneity of the scatter- 
ing cloud would further enhance the bremsstrahlung emission by a 
factor of q, see equation (|9}- 

It is also of interest to consider the dynamical state of scat- 
tering cloud. It cannot be free-falling as it would have given a 
too large accretion rate. It may be partly pressure supported. The 
dominant contribution to the pressure, P, is found to be radia- 



tion, P ~ T^LI{AnR^c). Then, the hydrostatic equilibrium, PjR ~ 
GMxiiipric^isI R^ yields Mx/Mq ~ tqL/L'^, where Lg is the Edding- 
ton luminosity for IMq (which is higher for He than that for H by 
the factor of /Je - 2). This yields Mx of several Mq. This is, how- 
ever, only an order of magnitude determination given our idealized 
model of a uniform cloud. Furthermore, the cloud is likely to be in 
a state of sub-Keplerian rotation, and thus partly centrifugally sup- 
ported, which would increase the required Mx- Thus, our model is 
consistent with the presence of a black hole in Cyg X-3. 

On the other hand, we note that our determination of the size 
of the cloud rests on the cut-off in the power spectrum. If that cut- 
off has a different origin, the scattering cloud could have a much 
larger size, and thus be effectively cooled by the stellar photons. 
Then, the models presented in Section [2l which are not consistent 
with the derived small size, could still be the correct ones. We note, 
however, that those models yield the temperatures of the intrinsic 
hot plasma substantially lower (and thus appearing less realistic) 
than both those of the models in Section |4] and those in the hard 
state of black-hole binaries. 

Interestingly, our calculations yield almost the same intrin- 
sic luminosities, Lx, for the two hard spectral states considered. 
In general, we expect softer sta tes to be more luminous (e.g., 
iDone. Gierliiiski & Kubota 2007). However, hysteresis effects in 
black-hole binar ies often lead to differ ent states to occur at the same 
luminosity (e.g.. lZdziarski et alj2004l) . and there are indication tha t 
such an eflrect also occurs in Cyg X-3 iSmit & van der Klislll996h . 
On the other hand, we have modelled only the hard part of the X- 
ray spectra. Given the very strong absorption in this system, there 
may be an additional soft component, stronger in the softer of the 
two considered states. An indication for that is the difference in 
the fitted values of A'h, with the apparently lower value of it in the 
softer state being possibly due to the presence of an additional soft 
emission. 

Our obtained value of Lx - 1.2 x 10'^ erg s"* is ~ 50- 
70 per cent of the values obtained for the corresponding spectral 
states by Hj09. Based on a comparison of the colour-luminosity di- 
agram for Cyg X-3 with that for black-hole binaries, Hj09 favoured 
Mx - 30Mq. Our lower value of Lx would imply a mass of ~ 2/3 
of that, provided this argument holds. Also, the parameters of the 
intrinsic X-ray emitting plasma derived by Hj09 would need to be 
significantly revised due to the existence of a scattering cloud pos- 
tulated in our model. 

An important issue is whether our model is applicable also 
to soft states of Cyg X-3 (3, 4, 5 in the classification of SZM08). 
Since Compton scattering of photons at a given energy moves them 
to a wide range of energies, there is no unique intrinsic spectrum 
corresponding to a given observed spectrum for given parameters 
of the scattering cloud, kTo and tq. Furthermore, the parameters 
of the cloud may change between states, in particular kTg follows 
from energy balance and thus depends on the incident spectrum. 
Then, a possible approach is to take observed soft state spectra of 
black-hole binaries and see how Compton scattering modifies them. 
However, this is a highly complex problem, given the wide range 
of observed soft state spectra, which req uires a relatively large 
number of free parameters t o model (e.g.. |Gierliiiski et alj 1 19991 : 
IZdziarski & Gierlinskill2004l : lOone et al.ll2007h . Thus, we defer its 
systematic study to future work. Here, we just consider one spec- 
trum, which we model in a very simple way. We take the average 
spectrum 3 and the cloud parameters fitted to the spectrum 2. As 
a simple intrinsic spectrum, we use the bmc model of xspec, which 
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Figure 7. A preliminary comparison of the spectrum of the average state 
3 (soft) of Cyg X-3, shown by the error bars, with a model consisting of 
an intrinsic emission of a Wien spectrum and a power law (dots) passing 
through a scattering cloud with the same parameters as those found for the 
state 2, To = 7.0, kTo = 2.5 keV. In the intrinsic spectrum, a power law 
with the photon index of 2 is convolved with a Wien spectrum with the 
temperature of 0.5 keV. The solid curve shows the intrinsic spectrum after 
transmission through the cloud, and the crosses show the transmitted spec- 
trum after absorption by external matter with A'h = 7 x 10^^ cm"-. For 
simplicity, neither bremsstrahlung nor Compton reflection/Fe K emission 
are taken into account. 



is a convolution of a power law with a Wien spectrun0. We find 
that this model gives a relatively good description of the average 
spectrum 3, see Fig.|7] Since we use a phenomenological, i.e., non- 
physical description of the input spectrum, we do not consider here 
the issues of energy balance or bremsstrahlung emission (which we 
defer to a future study with a physical model, e.g., that of hybrid 
Comptonization, see, e.g.. Gierliriski et al,,1999) . Still, our results 
show the potential possibility of explaining also the soft spectra of 
Cyg X-3 by a model with transmission through a scattering cloud. 

A test of our model can be performed by studying RXTE data 
for the hard states binned with respect to the orbital phase. From the 
expected anisotropy of the scattering cloud, the optical depth along 
the line of sight should be lower around the superior conjunction 
(the compact object in front of the donor) than around the inferior 
one. This would result in the cut-off of both the energy and power 
spectra being higher for the former case than for the latter. 



6 CONCLUSIONS 

We have considered the X-ray spectra of Cyg X-3 in its hard spec- 
tral state. The observed spectra show a cut-ofF/break at an unusu- 
ally low energy. We find this can be explained if a standard X- 
ray source is surrounded by a Thomson-thick cold plasma, which 
down-scatters the intrinsic X-ray spectrum. The presence of such 
a plasma was earlier postulated to explain the energy-independent 
X-ray orbital modulation and the lack of high frequencies in the 
power spectra. 



^ We note that although bmc is intended to model bulk-motion Comptoniza- 
tion in a converging flow around a black hole, it does not due to not includ- 
ing a high energ y cut-off, wh ich is always at ^100 keV for that process 
iNiedzwiecki & Zdziarskill2006) . Still, bmc is a useful generic approxima- 
tion to Comptonization by non-thermal electrons. 



We use a Monte Carlo code to model the efi'ect of Compton 
scattering in the cold plasma on both spectra and variability. For 
the latter, we obtain Green's functions and their Fourier transforms. 
Our results show that Compton scattering affects the power spectra 
at a frequency ~ 4 times lower than that obtained from considering 
the average number of scattering in the cloud. 

Finally, we consider thermal equilibrium of the cloud in the 
field of the incident photons. The obtained model has the cloud 
temperature of kT^ ^ 3 keV, the optical depth of Tq - 7, and the 
cloud size of R ^ 2 x 10' cm. We find that for those parameters, 
thermal bremsstrahlung emission of the cloud is important. Then, 
the cooling by this mechanisms facilitates achieving the exact en- 
ergy equilibrium of the cloud. The model reproduces both the X- 
ray spectrum and the power spectrum in the hard state. The cloud 
is partly supported by the radiation pressure. 



ACKNOWLEDGMENTS 

We than C. Done for valuable discussions, K. Leszczyfiski 
for help with the Monte Carlo code, and the referee for in- 
spiring suggestions. This research has been supported in part 
by the Polish MNiSW grants NN203065933 and 362/1/N- 
INTEGRAL/2008/09/0, and the Polish Astroparticle Network 
621/E-78/BWSN-0068/2008. We acknowledge the use of data 
obtained through the HEASARC online service provided by 
NASA/GSFC. 



REFERENCES 

Amaud K. A., 1996, in Jacoby G., Barnes J., eds, ASP Conf. Ser. 
Vol. 101, Astronomical Data Analysis Software and Systems V, 
p. 17 

Axelsson M., Larsson S., Hjalmarsdotter L., 2009, MNRAS, 394, 
1544 

Berger M., van der Klis M., 1994, A&A, 292, 175 
Davidsen A., Ostriker J. P, 1974, ApJ, 189, 331 
Dickey J. M., 1983, ApJ, 273, L71 
Done C, Gierliiiski M., Kubota A., 2007, A&ARv, 15, 1 
Gierliriski M., 2000, PhD thesis, N. Copernicus Astronomical 
Center, Warsaw 

Gierliriski M., Zdziarski A. A., Done C, Johnson W. N., Ebisawa 

K., Ueda Y, Haardt F, Phlips B. F, 1997, MNRAS, 288, 958 
Gierliiiski M., Zdziarski A. A., Poutanen J., Coppi P. S., Ebisawa 

K., Johnson W. N., 1999, MNRAS, 309, 496 
Gorecki A., Wilczewski W., 1984, Acta Astr., 34, 141 
Hertz P, Joss P C, Rappaport S., 1978, ApJ, 224, 614 
Hjalmarsdotter L., Zdziarski A. A., Larsson S., Beckmann V., Mc- 

CoUough M., Hannikainen D. C, Vilhu O., 2008, MNRAS, 384, 

278 

Hjalmarsdotter L., Zdziarski A. A., Szostek A., Hannikainen 

D. C, 2009, MNRAS, 392, 251 (Hj09) 
Illarionov A., Kallman T, McCray R., Ross R., 1979, ApJ, 228, 

279 

Kellogg E., Baldwin J. R., Koch D., 1975, ApJ, 199, 299 

Kylafis N. D., Klimis G. S., 1987, ApJ, 323, 678 

Kylafis N. D., Phinney E. S., 1989, in Ogelman H., van den 

Heuvel E. P. J., eds.. Timing Neutron Stars. Kluwer, New York, 

p. 731 

Lightman A. P, Lamb D. Q., Rybicki G. B., 1981, ApJ, 248, 738 
Magdziarz P, Zdziarski A. A., 1995, MNRAS, 273, 837 



© 2009 RAS, MNRAS 000,[TH9] 



Miyakawa T., Yamaoka K., Homan J., Saito K., Dotani T, 

Yoshida A., Inoue H., 2008, PASJ, 60, 637 
Nakamura H., Matsuoka M., Kawai N., Yoshida A., Miyoshi S., 

Kitamoto S., Yamashita K., 1993, MNRAS, 261, 353 
Niedzwiecki A., Zdziarski A. A., 2006, MNRAS, 365, 606 
Owocki S. R, Gayley K. G., Shaviv N. J., 2004, ApJ, 616, 525 
Poutanen J., Svensson R., 1996, ApJ, 470, 249 
Poutanen J., Lipunova G., Fabrika S., Butkevich A. G., Abol- 

masov R, 2007, MNRAS, 377, 1 187 
Poutanen J., Zdziarski A. A., Ibragimov A., 2008, MNRAS, 389, 

1427 

Predehl P, Burwitz V., Paerels R, Trumper J., 2000, A&A, 357, 
L25 

Rajeev M. R., Chitnis V. R., Rao A. R., Singh K. P, 1994, ApJ, 
424, 376 

Rybicki G. B., Lightman A. L., 1979, Radiative Processes in As- 
trophysics. New York: Wiley 
Shapiro S. L., Lightman A. P, 1976, ApJ, 204, 555 
Smit J. M. , van der Klis M. , 1 996, A&A, 316,115 
Sunyaev R. A., Titarchuk L. G., 1980, A&A, 86, 121 
Szostek A., Zdziarski A. A., 2008, MNRAS, 386, 593 (SZ08) 
Szostek A., Zdziarski A. A., McCollough M. L., 2008, MNRAS, 

388, 1001 (SZM08) 
van Kerkwijk M. H., 1993, A&A, 276, L9 
van Kerkwijk M. H., et al., 1992, Nat, 355, 703 
van Kerkwijk M. H., Geballe T. R., King D. L., van der Klis M., 

van Paradijs J., 1996, A&A, 314, 521 
Vilhu O., et al., 2003, A&A, 411, L405 

Vilhu O., Hakala P., Hjalmarsdotter L., Hannikainen D., Paizis A., 

McCollough M., 2007, in Grebenev S, Sunyaev R., Winkler C., 

The Obscured Universe. ESA SP-622, p. 377 
Vilhu O., Hakala P, Hannikainen D. C., McCollough M., Koljo- 

nen K., 2009, A&A, 501, 679 
Wardzinski G., Zdziarski A. A., Gierlinski M., Grove J. E., Jahoda 

K., Johnson W. N., 2002, MNRAS, 337, 829 
White N. E., Holt S. S., 1982, ApJ, 257, 318 
White T. R., Lightman A. P, Zdziarski A. A., 1988, ApJ, 331, 939 
Yamaoka K., Miyakawa T. G., Saito K., Uzawa M., Arai M., 

Sakamoto H., Yamazaki T., Yoshida A., 2006, in Proceedings of 

the VI Microquasar Workshop: Microquasars and Beyond. Proc. 

Sci.,p. 102 

Zdziarski A. A., Gierlinski M., 2004, Progr. Theor. Phys. SuppL, 
155, 99 

Zdziarski A. A., Johnson W. N., Magdziarz P, 1996, MNRAS, 
283, 193 

Zdziarski A. A., Poutanen J., Mikolajewska J., Gierlinski M., Ebi- 

sawa K., Johnson W. N., 1998, MNRAS, 301, 435 
Zdziarski A. A., Poutanen J., Johnson W. N., 2000, ApJ, 542, 703 
Zdziarski A. A., Lubihski P., Gilfanov M., Revnivtsev M., 2003, 

MNRAS, 342, 355 
Zdziarski A. A., Gierlinski M., Mikolajewska J., Wardzinski G., 

Smith D. M., Harmon B. A., Kitamoto S., 2004, MNRAS, 351, 

791 



Compton scattering in Cyg X-3 9 



© 2009 RAS, MNRAS OOO.dUl] 



